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, Abstract. We consider a modification of the coevolutionary voter model inspired by the idea 

of letting opinion dynamics reflect different strategies to promote consensus. This extension 
of the classic voter model breaks the symmetry of the two competing opinions, giving rise 
to more complex behaviour. A description of the phase diagram of the system is given, 
and the performance of the pair approximation is analyzed. The stochastic properties of the 
. active phase are studied, and the applicability of analytic approximations developed for the 

| coevolutionary voter model is discussed. 
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1. Introduction 

Many collective phenomena arise through structured interactions among a system's 
constituents HI HI- In this network of interactions, nodes (representing the system's agents) 
are assigned states and links to other nodes, with node states driving system dynamics along 
links. The coupling of node-state dynamics with their topological background is a widespread 
occurence and has been explored computationally as well as analytically in the field of 
adaptive networks [|3] 0] Q. 
^ ! The long-term behaviour is an important characteristic of adaptive networks and usually 

assessed by low-dimensional ODE models EH [3 [H 13. Possible asymptotic states are either 
1) frozen, for network configurations that prohibit further change or 2) active, if featuring 
perpetual dynamics in the thermodynamic limit. Frozen states in 1) are further classified 
according to how dynamics come to a halt: At absorbing consensus, a single node state 
is adopted globally and updating stops because of node-state homogeneity. The absorbing 
fragmented state on the other hand features status heterogeneity, but lacks active links that 
connect different nodes types and drive network dynamics. Adaptive networks in 2) can 
feature dynamic equilibria (DE), oscillatory behaviour or more complex dynamical regimes 
(see B21 SI for reviews). Moreover for finite-size networks in the active phase, stochastic 
fluctuations induce transitions from a prolonged active to a frozen state. The lifetime 
distributions of these metastable regimes can be described as first-passage problems, and the 
scaling of system size with several relevant averages can be uncovered fl6l ITOl ITTI . 

One simple example of an adaptive network with binary state space is the coevolutionary 
variant ifTOl of the classic voter model [12], mimicking the spreading of two opinions in a 
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population whose members strive for homophilic interactions. Its dynamics feed on active 
links connecting the two competing node types, and can be implemented in the following 
ways: In the node-update scheme, at every time step a node A and one of its neighbours 
B are randomly picked. If the link connecting them is active in the direct (reverse) node- 
update, then node A (node B) rewires it with probability cu to another randomly selected node 
of the same type, and adopts the state of node B (node A) with probability 1 — uj. If the 
link is however inert, i.e. connects two nodes of the same type, nothing happens. For the 
coevolutionary voter model with link update, a random link is picked. If it is inert, nothing 
happens. If it is active, one of its two end nodes is randomly chosen. Then that node rewires 
the active link with probability w to an arbitrarily picked node of the same type or adopts the 
state of the node at the other end of the link with probability 1 — cu. The ensuing dynamics 
are completely symmetric in node states, with the update schemes yielding similar outcomes 
[fT3l . yet with one major difference: For link update, the average network magnetization - the 
difference in the two fractions of nodes types - is conserved. Under node update however, it 
is just conserved if the network's degree distribution is sufficiently homogeneous, while it is 
not otherwise lfl4l . 

In the following, the symmetry of the coevolutionary voter model (referred to as 
symmetric model) will be broken by assigning different update schemes to different node 
ensembles and moreover by introducing a bias for the rewiring rule. The idea is that the two 
opposing opinions are associated with different social attitudes, which translates into different 
strategies to promote consensus in their holders' local environment. 

2. The Model 

In a network with state space {S, 1} and a fraction of I-nodes denoted by x, three elementary 
processes shall give rise to asymmetric adaptive opinion dynamics. 

In the first process, homophily is promoted through the asymmetric rewiring rule 
introduced in 0: S-nodes retract links from I-nodes with rate cu and reattach them to 
randomly selected S-nodes. Secondly, S-nodes shall adopt opinions from adjacent I-nodes 
with the transmission rate (1 — u)p (cu, p E [0, 1] without loss of generality). Both processes 
occur along randomly selected active links connecting S- with I-nodes, so that the S -ensemble 
engages in the coevolutionary voter dynamics of [10], but under a link-update scheme. 

The third process in contrast describes opinion adoption in the I-ensemble and is based 
on the node-update scheme: In guided relaxation, randomly selected I-nodes shall relax to S- 
nodes with rate (1 — co)(l — p)(l + am). The factor (1 + am) lets the network magnetization 
m = 1 — 2x E [—1, 1] steer the rate of the process, with a E [0, 1] quantifying the coupling of 
relaxation to m. Hence the relaxation of I-nodes is boosted by a strong (global) presence 
of S-nodes, whereas it is diminished by a dominance of I-nodes. Consequently opinion 
adoption among I-nodes resembles classic voter dynamics with (direct) node update, with the 
difference being that it is guided by overall network magnetization, not the magnetization 
of the neighbourhood of the respective I-node. It can therefore be seen as a mean-field 
description of the latter process. 
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These three processes set up an interplay of link-update coevolutionary voter dynamics 
with mean-field node-update voter dynamics by tying each to a specific node ensemble and 
letting those compete. The ensemble- specific update schemes shall reflect two competing 
spreading strategies in a population: Segregationist S-nodes are orthodox opinion holders 
that spread their opinion via social pressure and strive for local consensus by seeking to 
interact with their peers. Proselytic I-nodes engage with and convert S-nodes in personal 
interactions, and their heterodox opinion relaxes to the S-ground state at a rate that reflects 
the overall dominance of that opinion. While biased voter dynamics have been put forward 
through different rates of opinion adoption [15J or interactions along directed links [[T6l [TTIl . 
the asymmetric dynamics outlined above are a result of differing adoption strategies, their 
differing rates, as well as the biased rewiring rule described above. The dynamics are 
moverover reminiscent of the adaptive SIS model proposed in 0, in that in the latter I-nodes 
relax with a constant rate (retrieved through setting a = 0), while transmission and rewiring 
are as described above. 

A coarse-grained description of the network process is achieved by its moment 
expansion, where a set of evolution equations is compiled for the per-capita density of 
network motifs. In the subclass of pairwise models, only node and link densities are tracked, 
while at the same time higher-order motif densities essential to describe their evolution are 
approximated. This latter moment closure caps the hierarchy of motif equations at the level 
of links, featuring a parameter k that encapsulates the variance of the underlying network's 
degree distribution (see [fi~3l for a review). 

For the asymmetric opinion dynamics introduced above, the relevant link densities y and 
z refer to links among I-nodes (II-links) and to links connecting S- with I-nodes (STlinks), 
respectively. Since rewiring ensures link-number conservation, the density of links connecting 
S-nodes (SS-links) is given by (k)/2 — y — z for < y + z < (k}/2 and fixed mean degree 
(k). The standard pair approximation (PA) of the process, in the spirit of 0, then yields 



^ = (1 - u) {pz - (1 - p) r a x) 

|=(l-«)(^(«i4i + l)-2(l-rtr* 

dz z^ 
— = - z {u + (1 - u) (p + (1 - p) r a )) - (1 - u) pK- 



dt ' 1 — x 

+ 2(l-u)(l-p)r a y + 2(l-u) P K m ~ y ~ Z)Z (1) 

I — x 

with r a = 1 + a (1 — 2x). The last term of the time evolution of z in equations Q] for instance 
describes the gain in STlinks through the infection of either end of an SS-link. In that case the 
crucial density of the relevant triplet motif (consisting of a central S-node connected to both a 
S- and an I-node) is approximated using the densities of the aforementioned link types and of 
S-nodes. 

In the asymmetric opinion dynamics introduced above, link aquisition and opinion 
adoption are tailored to specific node ensembles, so that highly-skewed degree distributions 
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can ensue for a wide range of parameters (see []7] ED as examples for related dynamics). 
For that reason, instead of the regular random graphs taken in |[6l [TOj, we decide for initial 
Erdos-Renyi (ER) graphs (featuring a wider Poissonian degree distribution) and keep k = 1 
throughout, for it has been shown that for even more heterogenous degree distributions, the 
respective moment closure maintains validity [0 [T3). As we shall see, the PA sufficiently 
captures the corresponding network process with that choice of kappa. 

Choosing < a < 1, i.e. an intermediate coupling of guided relaxation to network 
magnetization, the same system phases as in the adaptive SIS model ensue (figure [T|). 
However the simple endemic phase shrinks considerably, while the simple bistable and stable 
oscillatory regimes are widened significantly. 
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Figure 1: Phase diagram of a) adaptive SIS b) asymmetric coevolutionary opinion dynamics with a = 0.99. Transcritical bifurcations 
(red line), fold bifurcations (green line), Andronov-Hopf bifurcations (blue line) and cycle-fold bifurcations (purple line) mark the following 
phases: S-consensus (S), simple active phase (A), simple bistable phase (AS), oscillatory bistable phase (OS, insets). Mean degree (k) = 5, 
compiled with [19]. 

Since within < a < 1 system dynamics is qualitatively the same as in the adaptive SIS 
model, we set a = 1 for all further discussions to fully couple the relaxation of I-individuals to 
the composition of the population, and refer to this case as the asymmetric model. Moreover, 
we will focus on the long-term behaviour of the dynamics and compare it to the asymptotic 
scenarios of the symmetric model. 

3. Asymptotic States in the PA 

Denoting the state vector of equations Q] as (x,y, z) with x E [0, 1], y E [0, (k)/2) and 
z E [0,(k)/2 - y] yields the two S- and I-consensus states (0,0,0) and (1, (k}/2, 0), 
respectively. It is straightforward to check that these are equilibria for equations CD A linear 
stability analysis of the PA reveals that the S-consensus is stable for 



2 -u 
~ (2 + <*»(l- W )- 
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As the PA equations are singular at (1, (k)/2, 0), one needs to resort to regularization 
techniques to obtain 

2 2 -3a/ 



p > mm 



! + (&}' 1-w 



as the parameter region for stable I-consensus (Appendix A i. Note that if < a < 1, I- 
consensus is not an absorbing state. 

Apart from a frozen phase, the PA yields a dynamic equilibrium (DE) for a small 
parameter region 

2-uj 2 

< P < — ttt (2) 



(2 + (k))(l-u) ~ r - l + (k) 
bordering the two consensus states. In this region of intermediate p and small u, relaxation 
and transmission balance out, with the small rewiring rate allowing for the continued existence 
of a nonzero density of active links and a steady-state fraction of I-nodes of 

2-(2+{k))p(l-u)-u 

UJ — p{± — LO) 

The system's only active phase ensues, whose size in the phase diagram shown in figure |2k 
decreases for increasing mean degree (k). For xa = 1/2, the magnetization of the system in 
steady state is zero, so that I-nodes recover with constant rate r and the PA describes the DE 
of the adaptive SIS model at the respective parameters. 
The triple point T in figure |2k at 

2 

= Pt 



l + (k) 
2 



uj = — - = w T 

3 + (k) 

lies at the confluence of the active and the two absorbing phases. It marks the end of the active 
phase and the advent of a bistable regime of coexisting consensus states for faster rewiring, 
with p tuning the competition of the latter' s basins of attraction. 



4. Comparison to the Stochastic Network Process 

Setting Initial and Exit Conditions 

To properly compare the PA with MC simulations, one has to faithfully translate network 
configurations into PA state vectors. An ER graph with a fraction x of randomly primed 
I-nodes is captured by the PA as (x , (k)xl/2, (k)(l — x )x ). Initial conditions in all MC 
runs are set this way and translated into the PA formalism accordingly, with the simulations 
employing the Gillespie algorithm ll2~0l . To let the PA capture the finite size of networks with 
iV nodes and N(k) /2 links, integration of equations |T] is stopped as soon as x < 1/N or 
1 — 1/N < x, emulating the two possible frozen consensus states in a network. 

Moreover it is important to identify "pathological" network configurations and avoid 
them altogether. If for instance MC simulations without rewiring ran on networks with 
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Figure 2: a) Phase diagram with steady states S-consensus (S), I-consensus (I) and the active phase (A), as well as the triple point T, 
bounded by solid lines. A sequence of symbols indicates the coexistence of respective attractors. The dashed line x\ = 1/2 yields DEs 
of adaptive SIS for respective parameter values, b) Change of asymptotic behaviour in the PA with initial conditions (0.1, 0.025, 0.45), 
(0.5, 0.625, 1.25) and (0.9, 1.025, 0.45) (numerical integration of PA, blue triangles) and initially connected ER graphs with fractions 0.1, 
0.5 and 0.9 of randomly assigned I-states (MC simulations, red squares). Regions of equal asymptotic behaviour are marked by the same 
sequence of symbols as in a). For sufficiently long simulation runs, stochastic fluctuations in the metastable state drive the full system into 
S- or I-consensus (inset of b)). Mean degree (fc) = 5, MC simulations with A' = 5000 nodes and results averaged over 100 realizations. 



isolated S-subgraphs (consisting of only S-nodes), I-consensus could not be reached, as these 
subgraphs would be left untouched by network dynamics. This would moreover add a constant 
offset to the network's magnetization that distorted relaxation. Consequently, "natural" 
system dynamics without rewiring can only unfold when taking place on a connected network. 
For consistency reasons, we ought to implement procedures that avoid fragmentation of 
initial random networks for all to, without ruling out later fragmentation in the course of 
coevolutionary dynamics. 

As (k) < log(iV) for mean degrees (k) and system sizes N used in MC simulations here, 
an initial ER graph is almost surely fragmented [j2T1 . Its isolated subgraphs need to be linked 
through a) linking two randomly selected nodes from separate subgraphs b) randomly picking 
a node that emanates links of the same type added in a), randomly choosing and deleting one 
of them c) repeating a)-b) until the graph is connected. Because the initial ER graph has a 
Poissonian degree distribution, a lower limit on its number of disconnected components can 
be given through e~^N, where e~( k ' is the average fraction of isolated nodes. For the values 
of (k) and N used, this lower limit approximates the actual number of initially disconnected 
subgraphs very well. Thus the fraction of nodes involved in this linking procedure is very 
small. The procedure neither introduces correlations in degree nor in status, so that apart from 
the vanishing isolated nodes, the main characteristics of an ER network are preserved. All 
initial networks used in the following MC simulations are connected this way. 

Comparing Asymptotic behaviour 

The possible coexistence of asymptotic states necessitates taking into account their competing 
basins of attraction. Selecting a sufficiently large set of initial conditions and monitoring 
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the resulting asymptotic behaviour of the system allows for the detection of all basins of 
attraction, both in integration of equations [Hand MC simulations. Browsing parameter space 
with this procedure would lead to a comparison of the phase boundaries of figure |2k with 
their MC analogue. Instead, we identify regions in parameter space for which a given small 
set of initial conditions lets dynamics drive the PA and the full system into the same set of 
asymptotic states. Comparing such parameter regions resulting from integration of equations 
[T]to those obtained from MC simulations allows for a quantitative comparison of PA dynamics 
and the corresponding network process without having to verify phase boundaries of figure 
Hi. 

This coarse-grained browsing of initial conditions over the whole parameter space yields 
a good agreement between numerical integration of the PA and MC simulations (figure 12b), 
indicating that for initial (connected) ER graphs, the PA faithfully models the actual dynamics. 
A S7-parameter region in that context means that either S- or I-consensus can be reached from 
the set of initial conditions used, whereas regions S, I and A signal a uniform asymptotic 
behaviour leading to S-consensus, I-consensus and a DE, respectively. In MC simulations, 
a metastable DE is observed for parameter values of the PA's active phase. In it, stochastic 
fluctuations eventually drive the system into one of the two consensus states. A more thorough 
PA description of the active phase, as well as a stochastic modeling of the corresponding 
metastable DE in the network process, will be given in section [5] and [6] 



5. The Active Phase in the PA 



For equations [T]in the active phase, the transient dynamics are reminiscent of what is reported 
from the symmetric model iflOl . in that the deterministic system relaxes quickly to a parabola- 
shaped slow manifold M D (figure |3k). In our case however, M D is generally not a line of 
equilibria, but spanned by two heteroclinics connecting the stable node (the DE) with the two 
saddles that represent the unstable consensus states. Once driven to M D , the system slowly 
advances along the respective heteroclinic towards the DE. 

The slow manifold Md is well approximated by a curve Ma given by 
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with Va{x} an d z^{x} being the equilibrium values of link densities y and z for an I-node 
fraction settling down to x. It follows that Ma is the set of all DEs that, for fixed (k) and 
cu < lut, are generated by all p for which the system is in the active phase (given by the 
interval in equation |2]). 

As p enters this interval from smaller values that lead to S-consensus, a transcritical 
bifurcation turns the stable node in (0, 0, 0) into a saddle, emanating a stable DE that is moving 
along Ma. It reaches (1, (k)/2, 0) at p = px an d, in another transcritical bifurcation, vanishes 
while turning the saddle there into a stable node representing I-consensus. 
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Since the vector field of the PA along M A is generally not tangent to the latter, M A 
usually does not exactly coincide with M D and is hence not a trajectory of the system. It is 
however straightforward to show that this matching improves for increasing (k) and tu, while 
it is already very good for the low mean degrees and rewiring rates considered here. As a first 
approximation, the description of the latter stages of system evolution towards the DE can, 
as in lfT0l l6ll. consequently be collapsed to one variable x, constraining the remaining two to 
yield a path along M A that does not depend on p. 

At the triple point T, i.e. for the highest rewiring rate still allowed in the active phase, 
the range of p for which there is an active phase shrinks to the single value p = Pt (equation 
|2]and figure |2k). It can be shown that then M A and Md exactly coincide, with the PA yielding 
a continuum of transversally stable stationary states given by equation HI A more thorough 
description of the model phenomenology at T will be given in section 17711 

6. The Active Phase in the Full System 

6.1. Identifying the Slow Manifold 

The onset of the metastable DE in MC simulations is characterized by vastly increasing 
convergence times needed to reach a consensus state, particularly for zero rewiring (figure |3j)). 
Projecting the network's time evolution onto the reduced phase space spanned by (x, y, z), the 
resulting random walk (RW) of the stochastic system follows closely the trajectory of its PA 
description along a curve M s (figure Hk), but with three fundamental differences: 

i) The exact shape of Ms depends on the system size N and approaches Md (the slow 
manifold of the PA) as N increases. 

ii) For various starting conditions well apart from M s , the RW does not immediately 
relax to Ms. Instead, the network undergoes a sequence of distinct configurations (corners of 
green trajectory in figure |3^) before its reduced description realigns with that given by M A . 
These intermediate network configurations are dependent on initial conditions and partition 
the RW into several segments. 

iii) Stochastic fluctuations drive the system along M s from the DE towards one of the 
two consensus states. Unlike in the stochastic process in the symmetric model, the type of this 
final consensus state in network dynamics is, for sufficiently large system sizes, predetermined 
by the system's parameters and independent of initial conditions (see below). 

6.2. Stochastic Route to Consensus 

Following the approach successfully carried out for the classic and coevolutionary voter model 
(6l [TOj, we will explore the possibility of reducing the description of the full system to a 
RW along Ms, while considering the simplest case of zero rewiring. Lacking an analytic 
expression both for M s and M D , the RW is assumed to take place along the approximate 
slow manifold M A given by equation 01 yielding a one-step process in the total number X of 
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Figure 3: The active phase in the asymmetric model, a) Time evolution of state variables x, y and z along M\ (black line) for u) = 0.05 and 
p = 0.32. All MC trajectories end up in I-consensus, starting from (0.01, 0.00025, 0.0495) (initial ER graph, red line) and (0.8, 0.2, 0.2) 
(maximally random graph with respect to initial conditions, green line). The latter initial network is obtained through i) generating two 
separate random S- and I-subgraphs compatible with given x, y and (k)/2 — y — z and ii) connecting them through N ■ z randomly 
assigned active links. Numerical integration from (0.8, 0.2, 0.2) (blue line) ends up in DE (x\, j/a, %a)- Network size N = 10 4 in MC 
run. b) Color-coded convergence times r in MC simulations within DE phase boundaries obtained from PA (solid green lines); maximum 
r are expected at x\ = 1/2 (dashed green line). MC simulations with N = 5000, averaged over 100 runs, c) System-size dependent 
splitting probabilities for u = computed through equation[6](triangles) and fixation probabilities taken from MC simulations (squares). 7Tt 
computed for p = 0.305 (starting from xq = 0.9, blue symbols), 7rg computed for p = 0.315 (starting from x = 0.1, red symbols), d) 
Convergence times in MC simulations as a function of system size for u: = and p = 0.3. Inset u> = 0.05. Results averaged over 10 4 runs 
(10 2 < N < 10 3 ) and 10 3 runs (10 3 < N < 10 4 ). Mean degree (k) = 5 in all figures. MC simulations in b)-d) from initially connected 
ER graphs, in b) and d) with xq = 0.5. Error bars in c) and d) are smaller than markers. 



I-nodes with the master equation 

= P Z A {X - 1}[X - 1] + 2(1 - p)(N - (X + 1))^V + 1] 

- P Z A {X} - 2(1 -p)(N- X)£[X\. (5) 

Here [X] is the time-dependent probability for a network of size N to have X = Nx I-nodes, 
and Z A {X} = Nz A {X/N} is the total number of active links computed from the respective 
link density at to = in equation 0] 

The consensus states X = and X = N are absorbing boundaries of the RW with 
starting point X e [0, TV], so that the eventual termination of the metastable state in the full 
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system corresponds to a first-passage scenario. One is interested in the splitting probability 
71-rpfo, N} (7r s {X , iV}) for the RW to end up at X = N (X = 0), approximating the fixation 
probability for I-consensus (S-consensus) in the full system. With equation \5\ and following 
Il22l . we calculate 



and 7r s {X , iV} = 1 — ni{X , N} for dynamics with u — 0. For a line of equilibria where 
transmission events balance out relaxation events, i.e. for pZ^{X} = 2(1 — p)(l — X/N)X 
at all X, equation [6] yields 7Ti{X , iV} = X /N as in the symmetric model. 

With equation [6l one can verify the earlier observation in MC simulations that for fixed 
X and increasing system sizes, fixation probabilities in the metastable regime converge to 
(for low p) or 1 (for larger transmission rates), indicating one intermediate p- value that would 
yield non-integer values throughout. Furthermore equation[6]reveals that for sufficiently large 
N, the final consensus type does not depend on X , but just on the position xa{w, p, n} of 
the DE on the slow manifold. This is plausible considering that for sufficiently large N, the 
system can be seen as initially drifting towards the DE, and from there diffusing to consensus. 
That also applies to to > 0, so that en route to the thermodynamic limit, stochastic fluctuations 
generally drive the system from the metastable to a consensus state whose type is increasingly 
independent of initial conditions. 

The increasing predetermination of the consensus type for large N is depicted in figure 
[3fc for both the full system and its reduced description along Ma. Even choosing initial 
conditions close to the disfavoured consensus state, its fixation probability, for large enough 
N, is observed to be proportional to e~ c N (with a constant c > 0), a scaling corroborated by 
equation [6l 

The preferred consensus for large N shall be determined by whether its distance to the 
DE, to be bridged in X by stochastic fluctuations, is smaller than N/2; an assumption that 
can be confirmed by equation |6] For xa = 1/2 then, i.e. when the PA describes the adaptive 
SIS model in steady state, one would conversely expect maximum first passage times for the 
RW along the slow manifold, and therefore also maximum convergence times r for the full 
system. Indeed this assumption yields a good estimate for maximum r for both the non- 
rewiring case and dynamics with topology change. The line of maximum convergence times 
in figure |3h differs slightly from x\(u, p, (k)) = 1/2 because the PA does not give the exact 
DE coordinates of the full system. Using in contrast the X-coordinates of the DE in MC 
simulations yields a very good match. 

Except for the qualitative considerations above, using equation[6]to compute tti{X , N} 
only yields moderately accurate predictions for the full system without rewiring (figure |3t), 
even when using Mg (as sampled from an ensemble of trajectories in corresponding MC 
simulations) instead of the approximate Ma. The reason is that since convergence times 
diverge quickly on static networks (figure |3]i), the computation of 7ri{X , N} is only feasible 
for relatively small system sizes. Yet for small N, assuming that the stochastic dynamics take 




-i 



(6) 
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place along a smooth M s neglects the effect of significant transversal fluctuations: Computing 
7Ti{X , N} along an averaged M s (taken over many stochastic trajectories) generally yields 
different results than directly averaging 7Ti{X , N} over the ensemble of trajectories. This 
remains an issue even for increasing system sizes, as M s retains "problematic" segments of 
low x and or z in the vicinity of the two consensus states, where transversal fluctuations are 
large even for large N. 

For symmetric dynamics along a line of equilibria as in [0 [KB, aforementioned 
limitations do not apply, since there fluctuations transversal to the random walk in X are 
decoupled from the latter. In our asymmetric model however, a stable DE is present on Mg 
with both the drift and diffusion of the random walk in X depending on the link density z. 
Even if Mg were a line of equilibria (section TT2\ . transmission and relaxation in its vicinity 
would generally not balance out, so that the (finite) full system with transversal fluctuations is 
not captured by its reduced description along Mg. Consequently a new framework is needed 
to address metastability in the asymmetric model. 

In the spirit of [|22| and similar to equation |6j one can derive an expression for mean 
first passage times for the process in equation |5] serving as an approximation for observed 
convergence times r in MC simulations. As however was demonstrated earlier, the system 
description generally cannot be collapsed to a one-dimensional RW in X, invalidating the 
derivation of further quantities. 

For to > 0, the sublinear scaling of convergence times r with system size N differs from 
that for co = (figure |3]1). Moreover it poses a stark contrast to the linear scaling found 
for the symmetric model [[TOl . considering that a stable DE on the slow manifold would be 
expected to counter diffusion towards consensus compared to dynamics on a line of equilibria. 
As asymmetric rewiring however increases overall degree heterogeneity, fluctuations in the 
aforementioned RW in X are enhanced, generally resulting in earlier consensus than on static 
networks or compared to dynamics with a symmetric rewriring rule. 

7. Comparison to Symmetric Coevolutionary Voter Dynamics 

Due to the symmetry imposed in the coevolutionary voter model of [[TOl . the network 
magnetization m is conserved in the thermodynamic limit for sufficiently homogeneous initial 
graphs [|T4l . and a single ODE suffices to describe the full system: It yields a (frozen) 
fragmented phase and an active phase which features a continuum of stable steady states. 
The specific steady state to be reached is then determined by the parameter m. For finite 
system sizes and due to stochastic fluctuations, the corresponding metastable state decays into 
either consensus state. It does so along a slow manifold given by aforementioned continuum 
of steady states. Consequently tti{X , N} = X /N, i.e. the fixation probabilities are 
system-size independent and scale linearly with starting coordinates X on the slow manifold. 
Therefore consensus in the symmetric model is reached from a metastable DE in the active 
phase - i.e. only stochastically, not dynamically. The model's symmetry encapsulated in 
the ODE moreover precludes any coexistence of phases; however stochastic bistability of 
the consensus states in the metastable state is given by the system-size invariant fixation 
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probabilities above. 

In contrast, the asymmetric model proposed here breaks the symmetry of the two 
node ensembles in a twofold way: Through state-dependent imitation rules and through 
tying the rewiring rule to just one specific node ensemble. The broken symmetry demands 
for additional degrees of freedom to capture the full system, yielding equations [Q As a 
consequence, network magnetization is not conserved anymore, but a system variable that 
guides relaxation. This in turn impedes a frozen state with < x < 1, so that the asymmetric 
model lacks a fragmented phase. The slow manifold is generally not a line of equilibria 
anymore, for isolated equilibrium points replace the continuum of steady states - one in each 
consensus state and, in the active phase, an additional stable DE. It follows that a) the steady 
state in the active phase is independent of initial conditions, particularly m b) consensus can 
now be reached dynamically c) consensus bistability is predicted by the PA and observed for 
the full system d) for consensus reached stochastically in the metastable state, the fixation 
probabilities are now system-size dependent, converging for large system sizes to either zero 
or one regardless of initial conditions. 

Unlike the symmetric model in |[T0l . the asymmetric model is not robust against varying 
topological backgrounds. Since changing network topologies may shift the balance in opinion 
competition, the initial topology can be crucial for the outcome of the asymmetric dynamics 
[|23ll . The (also asymmetric) adaptive SIS model on the other hand generally allows for 
sufficiently enduring coevolutionary dynamics to wash out initial differences in network 
structure [1241 H8l . This is because its active phase does not feature a slow manifold along 
which stochastic fluctuations towards consensus could be facilitated. 

When choosing for instance k, = ( (k) — 1)/ (k) in equations [T]to assume an initial random 
regular graph [fl3l . the intersection of the two consensus boundaries in phase diagram figure 
|2k disappears, and with it the active phase. Instead, bistable consensus stretches down to 
to = 0, featuring a slow manifold much like the active phase of the model with k = 1 did (see 
section |5]). This time however it is spanned by heteroclinics connecting two stable nodes in 
the consensus states with a saddle as the unstable DE. The different model phenomenology 
for that choice of k is corroborated by MC simulations, emphasizing the importance of initial 
topology for asymmetric opinion dynamics on adaptive networks. 

7.1. The Triple Point in the PA 

Tuning the rates of the asymmetric model's three elementary processes strengthens or loosens 
its asymmetry. As indicated in section |5l there exists a choice of parameters where symmetry 
is restored: At the triple point T that borders all phases, the PA in steady- state is exactly 
captured by the Ma in equation 0] (section [5]), such that for all < x < 1 

• Ma is formed by equilibrium points, as opinion- adoption events balance out according 

to pza{x} = 2(1 — p)(l — x)x. 

• (k$) = (ki) = (k) for mean degrees (k s ) = (2((k}/2 — y — z) + z)/(l — x) and 
(k{) = (2y + z)/x of the S- and /-ensemble, respectively. 
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These two steady-state equalities are not found anywhere else in the active phase of the 
asymmetric model and moreover yield z^\\ — x} = za{x} as well as ?/a{1 ~ x} = 
(k)/2 — Ua{x} — za{x} as a noteworthy consequence, i.e. at T in the asymmetric model 
and generally in the symmetric model, the system stays in DE even if all node states are 
flipped. 

It is straightforward to show that the line of equilibria given by M A is the same as for the 
PA of the symmetric model with link- update and without rewiring (given the same k — 1). 
Therefore at T, i.e. for a nontrivial choice of parameters, the steady states of the asymmetric 
model replicate the active phase of (one particular flavour of) the classic voter model. 

An additional feature unique to the asymmetric model contrasts topology change with 
opinion spreading: At T, transmission along and rewiring of Si-links happen at the same rate 
(1 — u)p = u. Combined with the aforementioned balance of transmission and relaxation 
events, this implies the equipartition of processes in steady state at T: Rewiring, relaxation 
and transmission then each account for exactly one third of events. As soon as rewiring 
dominates through w > wt, topology change impedes a dynamic equilibrium (figure |2k). 

7.2. The Triple Point in the Full System 

We find that for the approximate coordinates (wj ~ 0.150, p* T ~ 0.325) instead of (wt = 
\i Pt = §)> the full system with (k) = 5 features a drift-free attracting curve M s along which, 
on average, pz{x} = 2(1 — p)(l — x)x and magnetization is conserved (figure Hk). Indeed 
MC simulations suggest that T* is unique and marks the end of the active phase, similar to 
the triple point T in the PA. 

For (lo,p) approaching T* along the curve of maximum convergence times in figure 
3b, stochastic trajectories drift towards x = 1/2 on the curve Ms before slowly diffusing to 
consensus with increasingly small drift velocities. At T* the drift velocity is effectively zero, 
and for higher values of u the stochastic trajectories drift towards, instead of away from, the 
consensus states. This behaviour agrees qualitatively with the description given by the PA for 
the thermodynamics limit, in that the whole Ma is formed by equilibria when the change of 
stability of the equilibrium point on Ma occurs. 

While the shape of M s at T* encodes the dynamical aspect of neutral stability (the global 
balance of transmission and relaxation events), its stochastic features should also emerge, 
e.g. iii{Xo, N} ps Xq/N for all initial X considered. As T* is an isolated point, a high 
volatility of steady-state characteristics with respect to small parameter deviations is observed. 
Nonetheless fixation probabilities for the approximate T* do not display the usual convergence 
towards or 1 values for increasing system sizes (figure @J>). This is another indication, at the 
level of the stochastic properties of the system, of the vicinity of a point where transmission 
and relaxation events balance out along M s . 

In contrast to T in the PA description, the steady- state subensemble mean degrees of the 
full system are not equal at T* (inset of figure Sk). This is to be expected, as neutral stability 
and degree equality are independent constraints on M s , with the former occuring only at T* 
for a given (k). Consequently, asymmetric voter dynamics in the full system generally cannot 
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fully emulate the classic voter model in DE, featuring only a "weak" triple point T* as a ghost 
remnant of T in the approximate pairwise description. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

Figure 4: Identifying the triple point in the full system with plots for (ujJ, p^,) (red symbols) and (u) = 0.05, p = 0.3) (blue symbols) 
a) dynamically through balance of events AE = pz — 2(1 — p)(l — x)x and of mean degrees AK = (kg) — (hi) (inset) for bursts 
of simulations from fractions xq = 0.2, 0.5 and 0.8 of randomly primed I-nodes b) stochastically through the fixation probability for I- 
consensus as a function of xo for N = 100 (squares) and N = 1000 (triangles). Simulations in a) recorded for N = 10 5 and 10 < t < 100 
and in b) averaged over 10 4 runs and with error bars smaller than markers. All simulations from initially connected ER graphs and with 
mean degree (k) = 5. 



8. Summary and outlook 

A network process of two asymmetrically competing opinions is introduced, interpolating 
between adaptive SIS dynamics and a revolutionary voter model with two absorbing states. 
Only for the latter limiting case is the asymptotic behaviour found to be qualitatively different 
from the adaptive SIS model, and thus focussed on for the remainder of this work. The 
symmetry of the ensuing model is broken through ensemble-specific flavours of opinion 
dynamics: coevolutionary link-update dynamics then compete with mean-field node-update 
dynamics. 

Unlike the symmetric coevolutionary voter model, said asymmetric model is sensitive 
with respect to initial network topology. Choosing one particular initial topology whose 
pairwise description captures the asymmetric dynamics consistently, the resulting phase 
diagrams of the two voter models differ considerably: In the asymmetric case, fragmentation 
is absent, and a dynamical route towards consensus as well as consensus bistability ensues. In 
contrast, consensus can only be reached stochastically in the symmetric model's active phase. 
While an active phase is also present in the asymmetric model's PA, system magnetization is 
generally not conserved. 

In the full system, the active phase manifests itself through increasing convergence times 
towards consensus. As in the symmetric model, this metastable regime can be qualitatively 
characterized by reducing system dynamics to a random walk along a one-dimensional slow 
manifold. The latter is a line of equilibria in the symmetric case, where dynamics are trivial 
enough to preserve magnetization even during the system's transient. The presence of a stable 



An Asymmetric Co evolutionary Voter Model 



15 



DE in the asymmetric case entangles transversal with tangential fluctuations and complicates 
a quantitative stochastic description of metastability. 

The PA of the asymmetric model yields a triple point in parameter space where the 
system in DE emulates the classic voter model with link-update dynamics. For the full system, 
a similar point is identified that partially replicates the phenomenology of the triple point in 
the PA: While it features an attracting curve along which magnetization is conserved, the two 
node ensembles in DE do not display equal mean degrees. However the adaptive SIS model, 
for chosen parameter combinations, yields node-ensemble symmetries in steady state that are 
also displayed by the full system: The respective PA features (k s ) = (hi) at (1 — to)p = u> as 
does our asymmetric model, while in contrast this equality not only holds in MC simulations, 
but extends to steady-state subensemble degree distributions [18J. 

In previous work it was shown that, given certain ergodic properties of a coevolutionary 
process in DE, steady-state averages imply steady-state distributions they arise from, and that 
these steady-state distributions can be obtained through an analytical framework [|2~5l IT8l . 
The latter could hence be used to characterize the steady-state subensemble networks in the 
metastable regime of the asymmetric model. 

Preliminary results moreover indicate that for the symmetric model, degree distributions 
remain unchanged along the line of equilibria. The aforementioned framework, initially 
designed to detect and describe pointlike DEs, could be extended to analytically support this 
observation. 

We furthermore leave for future work a satisfactory quantitative description of the 
stochastic properties of the metastable state. Encouraging analytic results in that direction 
have been obtained for fluctuations around the DE of a slow manifold [|26l . However to fully 
characterize metastability, one would need to extend existing approaches to the entire slow 
manifold. 
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Appendix A. I-Consensus as a Fixed Point 

Expressing equations \T\ in spherical coordinates 9,(p G [0,7r/2] and r e [0,1] with x = 
1 — rsin(0) cos(0), y = (k)/2 — rsin(6 l ) sin(0) and z = rcos(9), one is interested in the 
flow's behaviour in the vicinity of I-consensus, i.e. for r — > 0. 

In these coordinates, the singularity at I-consensus can be regularized, and we find that 
r = for r = and all angular coordinates. Hence the sphere octant E at r = is invariant 
and represents I-consensus, the latter of which is then indeed a steady state. Parameter 
regions of stable I-consensus can consequently be determined by identifying stable fixed 
points (9*, <p*) of the angular flow on E, in conjunction with demanding stability of the radial 
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flow at (#*, <p*). In the following, "fixed point" will be exclusively used in the context of the 
two-dimensional angular flow on E, whereas I-consensus in the full system is referred to as 
such. 

It is straightforward to classify the parameter dependence of the angular flow into six 
typical phase portraits, each of which yielding three fixed points at the boundary <\> = n/2 and 
none, one or two in the interior of E. With these phase portraits at hand, the stability of fixed 
points is separately investigated for those two regions. 

At the boundary <p — 7r /2, the angular flow has fixed points (#*,</>*) at 
arccos {2/v / 5},tt/2), (7r/4,7r/2) and (tt/2, tt/2). The first one is unstable, whereas the 
remaining two are non-hyperbolic with one zero and one negative eigenvalue each. The radial 
flow is unstable at (ir/2, tt/2) and stable at (7r/4, 7r/2), so that only the latter fixed point is of 
further interest. Abovementioned phase portraits reveal it to be stable as soon as there are two 
(for p < 2/(1 + (k))) or no (for p > 2/(1 + (k))) fixed points in the interior of E. For both 
cases, this criterion translates into a comparison of tangent slopes of respective nullclines at 
(tt/4, tt/2) and yields 

p> (2-3w)/(l-w) 

as a sufficient condition for stable I-consensus. 

Consequently for p < (2 — 3u) — there is exactly one fixed point in the interior of 
E, and it follows from aforementioned phase portraits that it is always stable. It can moreover 
be shown that the radial flow at this fixed point is stable if 

p> 2/(1 + (*;», 

so that overall, I-consensus is stable for 

p > min{2/(l + (k)), (2 - 3w)/(l - w)} . 

Note that there can be coexisting stable fixed points for the angular flow - one at 
(tt/4, 7r/2) (automatically implying stable I-consensus) and one in the interior of E. For 
the (k) chosen here, the basin of attraction of the former is significantly smaller than that of 
the latter. If the radial flow is unstable at the interior fixed point, trajectories in the vicinity of 
I-consensus usually get attracted to those angular coordinates and are consequently repelled 
from I-consensus in the r-direction. This is the case for low p and large to, so that reaching 
I-consensus for that parameter region requires careful selection of initial conditions (figure 
12b)- If however the radial flow is also stable at the interior stable fixed point, then there are 
two routes towards I-consensus, characterized by the two different angles under which the 
trajectory approaches it. 
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